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Abstract. The recent Sumi et al (2010, 2011) detection of free roaming planet 
mass MACHOs in cosmologically significant numbers recalls their original detection 
in quasar microlening studies (Schild 1996, Colley and Schild 2003). We consider 
the microlensing signature of such a population, and find that the nano-lensing 
(microlcnsing) would be well characterized by a statistical microlensing theory 
published previously by Refsdal and Stabel (1991). Comparison of the observed First 
Lens microlensing amplitudes with the theoretical prediction gives close agreement 
and a methodology for determining the slope of the mass function describing the 
population. Our provisional estimate of the power law exponent in an exponential 
approximation to this distribution is 2.98+J;g. where a Salpeter slope is 2.35. 
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1. Introduction 

The recent detection of planetary mass MACHOS seen in high-cadence searches toward 
the Galactic Center and the Large Magellanic Cloud (LMC) (Sumi et al, 2010, 2011) 
suggest that a significant population of dark planet-mass MACHOs populate the halo of 
our Galaxy, which may constitute the Galactic dark matter. This would (partly) explain 
the "missing baryon problem", the fact that about 90% of the baryons in the solar 
neighborhood are unaccounted for. Indeed, below we clarify this from similar detections 
made in last decades with micro-lensing and even nano-lensing in the Q0957+561 A,B 
gravitational lens system (the so-called First Lens). Insofar as the observed MACHOs 
are planetary mass bodies dominated by hydrogen and 26% in weight of helium, other 
modes of their detection as extreme scattering events seen in the line of sight to 
ordinary quasars, and as the sources of the ubiquitous "dust" emission at temperatures 
of minimally 15 K and of the "mysterious radio events", may also be signaling their 
detection (Nieuwenhuizen et al 2010). 

In the following we re-consider the expected microlensing signature, optical depths, 
and mass spectrum of the observed planet mass MACHO population. Presently there is 
confusion in the literature about what to call these objects, because there is increasing 
awareness from emerging statistics of ubiquity of orbiting planets accompanying 
ordinary stars that the many orbital interactions should frequently result in escape 
from the initial orbiting planetary systems. To avoid confusion, we call free-roaming 
planetary mass condensed objects "escapees" if originally formed in pre-stellar accretion 
discs, and "micro brown dwarfs" (/xBDs) if formed primordially by gas fragmentation. 
Of course we still call ordinary star-orbiting bodies planets. 

1.1. Quasar Microlensing and the missing baryons problem 

The existence of a population of planet mass microlensing objects, also called MACHOs, 
was first inferred from quasar microlensing studies by Schild (1996) when it was 
discovered that no value close to the accepted time delay would remove the pattern of 
daily sampled quasar brightness fluctuations. The subject was complicated by the fact 
that to securely recognize the signal, it would be necessary to measure the double image 
quasar time delay to a fraction of a day. This was accomplished by an international 
around-the-world consortium (Colley et al, 2002, 2003) whose 417.1-day time delay 
value still stands as the most precise time delay ever measured. Re-analysis of an earlier 
observation of 5 consecutive nights of continuous brightness monitoring data produced 
a microlensing event of duration 5 hours (quasar proper time), by Colley and Schild, 
(2003). 

Such rapid events can realistically be understood only as resulting from 
microlensing. Recall that a quasar is approximately 100 times brighter than our Milky 
Way galaxy. If attributed to the quasar, the observed 1 percent brightness fluctuation 
in the quasar light is energetically equivalent to the entire Milky Way luminosity being 
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switched on and off in only 5 hours - which is an implausible scenario. [|] For observation 
from the earth, the intensity fluctuation can be accomplished by the gravitational field 
of a microlensing compact object. Indeed, an observed negative brightness change is in 
fact a re-direction of the 1% part of the enormous luminosity away from the line of sight 
of the affected quasar image. Likewise, a positive 1% change is due to a lensing effect 
that focuses more light in the observer's direction. 

The original rapid microlensing detection was greeted with strong interest, because 
it was immediately recognized that it could not be caused by stars (Catalano 
1997). Since it was obviously seen as a continuous brightness fluctuation pattern at 
approximately unit optical depth, it must represent the detection of the dark baryonic 
matter. But the interest faded away when the EROS and MACHO consortia did not 
observe similar MACHOs in front of the Magellanic Clouds. However, the quasar data 
themselves have never been questioned and related effects were observed on other, lensed 
quasars, as commonly discussed in the context of measurement of time delay (Burud et 
al, 2000, 2002; Paraficz et al, 2006, Vakulik et al, 2007, 2008). To settle this dispute, 
we intend to redo a search in front of the Large Magellanic Cloud (Schild et al, 2012). 

To explain the quasar observations, other possibilities than MACHOs were 
considered as well, such as hypothetical orbiting luminous blobs in the accretion disc, 
for which, however, there had never previously been evidence (Gould and Miralda- 
Escude, 1996). Finally, a series of simulations of orbiting luminous blobs and obscuring 
clouds by Wyithe and Loeb (2002) produced brightness curves that could be compared 
to observations. Simulations for orbiting dark spots and bright spots microlensed by 
patterns of cusps originating in the lens galaxy show that the longer duration events 
have smaller brightness amplitudes than shorter duration events (Figure 6). But the 
wavelet analysis of the observed microlensing brightness fluctuations by Schild (1999) 
consistently showed an opposite pattern of larger brightness fluctuation amplitude for 
longer duration events. Moreover, in general the simulated brightness curves do not 
look like the observational results, and in particular do not show the observed feature 
of equal positive and negative events (Schild 1999). 

An additional simulation showed that for the process championed by Schild (1996), 
with quasar structure microlensed by cusps originating in the lens galaxy, larger 
amplitude brightness effects were always found for longer duration events (Wyithe and 
Loeb, 2002, Fig. 9). The wavelet analysis of the microlensing brightness fluctuations by 
Schild (1999) consistently showed this pattern of larger brightness amplitude for longer 
duration events. 

Since the Wyithe and Loeb (2002) simulation only covered stellar to Jupiter 

| There is a counter indication for this argument. Colley and Schild (2003) report a 5 night observation 
for both the A and the B image in periods corresponding to a common quasar time. In each day (night) 
of the 5 day observation by Colley and Schild, there is a common fluctuation in both images of the 
quasar with period of approximately one day. On the last night the trends in both images are different, 
which is interpreted as an additional lensing event on top of a common trend. However, the 24/7 
monitoring of the quasar for 10 days in 2001 did not exhibit this trend (Colley et al, 2003). This 
conundrum has to be resolved with new observations. 
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microlensing masses and not to the Earth masses or smaller, as implied by the Colley 
and Schild (2003) event, the targeted simulation by Schild and Vakulik (2003) for The 
First Lens, Q0957, is of more relevance. It demonstrates microlensing events caused 
not by orbiting blobs, but rather by a drift pattern of microlensing cusps originating 
in the lens galaxy and microlensing the discrete luminous quasar structures, creating a 
continuous pattern of brightness cusps at the observed 1% level with event durations of 
a few days, caused by microlens masses 1O _5 M . 

Thus in this paper we return to the interpretation of the quasar microlensing effect 
as due to a population of MACHOs. In the following sections we will first consider 
in section 2 the statistics describing the probability of microlensing in the quasar lens 
system, and then in section 3 we show how the theory of microlensing for large luminous 
sources gives an excellent fit to the observed amplitude of brightness fluctuations. In 
section 4 we demonstrate a technique for determination of the mass function describing 
the observed masses of the rogue planet MACHOs. We summarize our conclusions in 
section 5. 

2. Quasar structure and microlensing observations 

2.1. On size scales of luminous quasar structure inferred from microlensing 

Technically speaking, the subject of lensing by individual objects within a resolved 
object like a distribution of stars in a lens galaxy foreground to a distant quasar 
is discussed under the topic of microlensing, although if the grainey distribution is 
composed of planet mass objects, it would more correctly be described as nano-lensing. 
We follow standard usage and simply adopt the word microlensing in this report. 

It is also important to understand that the microlensing by planet mass objects 
resolves the structure of the central regions of quasars, taken to be black hole or MECO 
objects (Schild, Leiter, and Robertson 2006; SLR06). For one such object, the doubly 
imaged quasar Q0957+561 the distance observed between the A and B images is 6.26 
arcsec. It has redshift of source z$ = 1.43 and of the lens zl = 0.355, the BH mass is 
3.6 1O 9 M (SLR06) and the gravitational radius R G = GMc~ 2 is 5.3 • 10 14 cm. 

We describe the luminour ring at the inner edge of the accretion disc as a torus 
having two radii of importance; the outer radius, which refers to the quasar central 
radius to the inner edge of the accretion disc, and the sectional radius, which is half the 
thickness of the accretion disc and therefore the cross sectional radius. Then the outer 
radius is -Router = 74i?c = 3.94 10 16 cm and the sectional radius is R' ± = 1.52 • 10 13 
cm. The main part of the quasar light comes from a structure of this radius. [§] 

§ Indeed it is known from the examination of the quasar spectrum that at this wavelength of brightness 
monitoring, the inner edge of the accretion disk contributes only ~ of the observed brightness. What 
is microlensed is a bump in the spectrum, which is the thermal peak of light originating at the inner 
accretion disk. It is known that the brightness of this structure is \ of the total brightness at the 
monitoring wavelength. So if the | of brightness coming from the outer larger structure were absent, 
the quasar would be fainter overall, and the observed amount from the microlensed structure would be 
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We adopt the standard cosmology with H = 70 km/s Mpc, fi c = 0.23, £Ib = 0.05, 
so that Qm = 0.28, and Q\ = 0.72. The general formula for the angular distance 
oIa{z) = Hq 1 ^ + z)~ x Jl +Z xx~ 1 {VLl + ^m^ 3 ) ^ 2 yields for the angular distance to 
the source ds = cIa(zs) = 1160 Mpc and to the lens di = = 890 Mpc, 

while the angular distance between them, for light that we observe now, is d$L — 
ds — (1 + z L )(l + z s )~ l d L = 670 Mpc. The angular radius of the source is thus 
Pouter = R outer /d s = 2.14 //as. 

The Colley and Schild (2003) observation of an event of duration 12 hours 
(observer's clock) as an event with an approximately Gaussian profile instead of a 
more sharply peaked brightness profile (as from a point source) can be used to give an 
approximate dimension to the luminous quasar structure. Assuming as in SLR06 that 
the brightness profile results from a microlensing rogue planet moving past the luminous 
inner edge of its accretion disc, for a standard cosmological transverse velocity of 600 
km/sec, and a quasar distance of 1.31 times the lens distance (in the above mentioned 
cosmology) we conclude that the ring-shaped luminous inner edge of the accretion disc 
has a thickness dimension approximately equal to the Einstein ring diameter of the 
microlens, or d' ± = 1.52- 10 13 cm for the Colley and Schild (2003) event in the cosmology 
Hq = 75 km s _1 Mpc -1 , Vl M = 1, = (as indicated by the prime on R' ± ). Hence 
its angular diameter from that cosmology is 9± = d' ± /d' A {zi) = 1.41 nas, which in our 
cosmology corresponds to a radius d± = 1.88 10 13 cm. The black-hole-centric diameter 
of the inner luminous edge of the accretion disc has been determined from reverberation 
to be 7.8 • 10 16 cm (SLR06; Schild and Leiter 2009), so multiplying by 2n times R± 
determines the area of the luminous structure as 2.83 • 10 30 cm 2 . In the theory of RS91 
the number of deflectors is the product of their normalized surface density a and the 
surface area, so we may replace the luminous ring by an equivalent round structure 
having this area for an effective radius of R eS = 0.949 10 15 cm and effective angle 
9 cS = R cS /d L = 71 .4 nas. 

It is necessary to calculate this luminous area radius to compare with the radius 
of the Einstein ring of a microlensing particle. For the Colley and Schild (2003) 
microlensing event, the radius of the Einstein ring was previously calculated from 
the event duration to be | • 1.52 • 10 13 cm = 0.76 • 10 13 cm and therefore the RS91 
requirement that the luminous structure's radius should be at least 5 times greater 
than the microlens's Einstein radius is well satisfied. Therefore we may use the RS91 
statistical theory results to describe the relationship of brightness fluctuation amplitude 
to microlensing optical depth. 

This allows us to use two statistical results that relate the measured brightness 
fluctuation amplitude to the area of the lensing luminous source and to the mass function 
of the nano-lensing objects. 

the same. But the microlensing comes from j of the quasar brightness. Since we observe fluctuations 
in the total quasar brightness, the observed fluctuations should be multiplied by 4 to get the relative 
brightness fluctuations that would have been observed if the outer luminous structure were absent. 
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3. Estimation of the Einstein ring radius and microlensing event duration 

The discussion until now pre-supposes that all of the rogue planet MACHOs constituting 
the baryonic dark matter have the same mass. That mass was predicted to be typically 
1O~ 7 M = 0.03M e by Gibson (1996), who also predicted that the initial fragmentation 
would have been immediately followed by an accretional cascade to larger masses, 
especially in the times immediately following fragmentation after recombination, because 
the MACHOs were hotter and larger from their gravitational energy release, and because 
the density of the universe has been monotonically decreasing since the formation 
epoch, increasing the spacings. A more recent estimate of the mass comes out larger, 
3.9 1O" 5 M = 13M e (Nieuwenhuizen et al, 2009). 

Thus the mass function F(M) describing the number of MACHOs by F(M)M as 
a function of mass would have been modified since the original formation epoch by the 
process of accumulation. The prescient statistical microlensing theory of Refsdal and 
Stabell (1991) that describes the determination of brightness fluctuation amplitudes in 
terms of the Einstein radius of the MACHO and the optical depth, also demonstrates a 
formalism for determining the slope the mass function F(M). 

This occurs because relative to a simple mass function with a delta function for 
some mass M, a mass function with additional masses larger than M has relatively 
more large-amplitude events. We shall consider the case of F(M) expressed as a power 
law, 



while F = for M < Mi and for M > M2. In practice, the mean amplitude of brightness 
fluctuations is measured from brightness monitoring, the Einstein Ring diameter is 
estimated from the duration of the microlensing events, the diameter of the source is 
known or estimated, and the optical depth is known from the overall macro-lensing 
model. 

We adopt a standard transverse velocity of v = 600 km/s, based upon the 
presumption that extreme cosmological departures from the co-moving expansion 
velocity are top-limited by the local Great Attractor at approximately 1000 km/s, and 
that other velocities, (e.g., Earth orbit, Galaxy rotation, Galaxy cluster motion, etc) 
are of order 50-250 kms/s and uncorrelated. So a value considered accurate to a factor 
\/2 uncertainty, is ordinarily taken to be 600 km/s. 

The opening angle of the Einstein Ring is 



which yields 1.1 10 11 rad = 2.3 yuas for M = M & , explaining the name "microlensing" 
for solar mass objects and 4.0 nas, "nanolensing", for objects of earth mass. The M Q - 





(1) 




(2) 
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Einstein radius is Re = Qe&l = 3.05 ■ 10 16 cm, (Refsdal et al, 2000), the microlensing 
cusp crossing (proper) time is 16.1 year = 5880 days. 

From Eq. (J2| we obtain in Table 1 some fiducial values of microlens masses and 
microlensing event times. In particular we find that the observed 12 hr event detected 
by Colley and Schild (2003) results from a nano-lensing mass of 2.7 • 1O _9 M . 



mass 


duration (days) 


M 


5880 


0.01M© 


588 


M Jup 


182 


O.OOO1M 


58.8 


-^Earth 


10.2 


10- 6 M 


5.9 


10- 8 M 


0.59 


2.9 ■ 10~ 9 M Q 


0.2 



Table 1: event duration (proper time) as function of the lensing mass. 

The bottom line of Table 1 is for the shortest quasar microlensing event ever 
observed, (Colley and Schild, 2003) with an observed event duration of 12 hours and 
a cosmologically corrected event duration of 5 hours (quasar local clock). Because it 
cannot be considered certain that the observed rapid event is due only to a nano-lensing 
cusp crossing, but possibly due at least in part to some effect of orbiting luminous 
blobs or obscuring matter (Gould and Miralda-Escude, 1997), the case for Lunar- 
mass detection is not certain; however the observed brightness record for the 5 days 
preceding the securely detected event allow for the detection of more events of similar 
low brightness amplitude and durations of only hours (Colley and Schild, 2003, Fig. 1). 

4. Estimation of the rogue planet mass function 

A statistical microlensing theory for microlensing by masses with Einstein Ring 
diameters smaller than the size of the light emissing structures has been given by Refsdal 
and Stabell (1991). In its simplest implementation, it shows the expected average rms 
brightness fluctuation amplitude expected for a random distribution of microlenses all 
assumed to have the same mass (RS91, Eq. 1). A further elaboration of the theory 
also considers the case of microlensing by a distribution of masses having a power law 
distribution function as normally assumed for stars. In the latter case, the exponent is 
found to be a = 2.35 and the stellar mass distribution is called the Salpeter function. 

Thus if the optical depth is already known, as from a detailed model of the macro- 
lensing producing the double quasar image, then a correction to the measured brightness 
fluctuation mean amplitude can be determined and from the RS91 Fig. 2 plot, a 
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correction for the slope a of the mass function F(M) determined, if the ratio of the 
upper and lower mass bounds can be estimated. 

In our implementation of this scheme, we adopt the following parameters. We adopt 
an optical depth of a = 0.707 (Refsdal et al 2000). The area of the luminous quasar 
structure involved in microlensing was estimated in section 2 to be 3.67 • 10 30 cm 2 for 
an equivalent radius of 1.08 10 15 cm. It is necessary to calculate this just to ensure that 
the ratio of radii is larger than 5 for the RS91 statistical theory to be applicable, and 
we find a ratio of 6.7. Thus application of the RS91 theory is appropriate. 

We shall need integrals of powers of M, 



/ M F(M) M k = —. 2 1 = \ . 

hh ' Ml~ a k + l-a Ml~ a k + l-a 



with X = MijM\ > 1. The prefactor A can now be fixed by the normalization Iq = 1. 
The average of M k is 



k _ J^MF(M)M k _I k 
{M } ~ /JMF(M) " V (4) 

which does not depend on A anyhow. RS91 define the effective mass: 



M cS = ^ = MF 2 (a) (5) 
with M = y/M~iM~2 the geometric average of the upper and lower mass. Furthermore, 

/ J sl/2 i 

F{a) = {MTj =xvi 

RS91 point out that = 1 for all X. They then argue and support by simulations 
that for such a mass distribution, the predicted rms amplitude of the brightness 
fluctuations for microlenses of effective mass M e e reads 



2-a X 



3-q 



3 - a X 2 - a - 1 



(6) 



X oi7 1/2 E (M eS ) 

d)m R Q = 2.17 • cr 1 ■ 



0_ 

= 2.17 • a 1 / 2 ■ MM) . F ( a ) (7) 

u 

where 9 is the angular radius of the source. As discussed above, we shall take 
6 = 9 cS = i? e ff/^L = 71 nano-arcsec. 

The critical observational result determining the mass function exponent a is the 
observed rms amplitude of the observed brightness fluctuations (the value of 5m obs in 
equation rtTJ)). A plot of the rms brightness fluctuations for the Q0957 gravitational lens 
system revealed that a linear relationship exists between the event duration and the 
amplitude of the measured fluctuations. This is represented in Fig. 8 of Schild (1999) 
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and the best fit curve describing the rms amplitude as a function of event duration t ev 
was given as: 5m ohs = 0.00133 t ev day" 1 + 0.0060 mag. The measurements extended 
over a range of t ev -values from 2 to 64 days. However it was shown in Schild (1999) 
that the quasar has fine inner luminous structure presumed to be responsible for the 
fluctuations discussed here, and a more diffusive outer luminous structure. However 
in SLR06 it was shown that the measured quasar brightness fluctuations were at an 
observed wavelength of 680 nm which for cosmological redshift 1.43 originated at 280 
nm which is at the "small blue bump" in the quasar energy distribution. Moreover the 
fraction of quasar luminosity originating in the small blue bump is 1/4 of the total, so 
a correction factor 4 must be applied for the dilution of the microlensed inner structure 
by the outer UV-optical continuum measured in reverberation as described by Schild 
(2005) and also as modeled by Schild & Vakulik (2003). 

With correction for the outer quasar luminosity the inner microlensed region's rms 
brightness fluctuation amplitude is (see also the footnote at pages 4-5) 



We have determined our fitting interval with limits Mi and M 2 as follows. The 
lower mass limit inferred from the hydrodynamical theory for the formation of small 
microlensing particles predicts a lower mass limit of Mi = 2.7 1O _7 M which corresponds 
to a wavelet duration of 2 days on observer's clock, the lowest value for which we 
have a direct measurement of mean amplitude from detection of approximately 220 
microlensing events in the 4-year intensely sampled time interval analyzed by Schild 
(1999). The upper mass limit is estimated to be M2 = 2.4 • 1O~ 4 M from the maximum 
of the wavelets observed. 

For this range of observed masses, the effective mass reads M e g = (Mi • M 2 ) 1 ^ 2 = 
8.05 ■ 1O" 6 M = 2.7M and X = M 2 /Mi = 890. For a microlens of this mass the event 
duration would be 3.67 days from the calculation of its Einstein ring diameter and for 
the assumed microlens transverse velocity as above. 

The event duration that enters (l8l is defined for each deflector mass as 



with v = 600 km/s. In our statistical consideration we should average this over the 



Sm ohs = 0.00532 t ev day" 1 + 0.0240 mag 



(8) 





mass distribution Q. This yields 




26 E (M)d L 



G(a) = 8.28G(a) day, 



(10) 



where 



G{a) 



a - 1 1 - X 3 / 2 



a 



(11) 



a-\ l-X l - a 
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Figure 1. Upper curve: H(a) for the case X = 890. The maximum 6.09697 occurs 
at a = 1.85815; the minimal value H — 1 is attained for a — > ±oo. Lower curve: the 
right hand side of Eq. (12 1. The curves intersect at a* = 2.98. 



We can now equate {Sm ^ s ) to ^m^g, which can be written as 



H(a) 



H obs + H . 



(2) 
obs 



H, 



(i) 



obs 



0.0049(1 + z L )d L e eS 
y/av day 



G(a) 



1.42, H, 



(2) 
obs 



0.00276 OesX 1 ^ 



0.774 



where 



H{a) = X 1 / 4 



F(a) 
G(a) 



a-l 1 



X 



l-a 



a 



1 1 - X 3 / 2 - c 



a 



1 -X 3 



3 - a 1 -X 



2-a 



1/2 



(12) 



(13) 



(14) 



For our value X = 890, H has a maximum Hr, 



6.097 at a r = 1.86, while H -> 1 for 



a — > ±oo. A plot of both sides of Eq. (12) is presented in Figure 1. They intersect at 
a* = 2.98 and we estimate the error as a* = 2.98 +10 — 0.5 . We neglect the intersection 
for small a, so that we have in any case a > a c = 1.86 and, in particular, we may end 
up close to the Salpeter value a = 2.35. However, there is no reason to expect the two 
values to agree, since they are likely to be dominated by different processes in star and 
in primordial rogue planet formation. 



We list in Table 2 the numerical values determined for the constants in Eq. (12). 



The value for the microlensing optical depth a is the image A value from Refsdal et al 
(2000) from the image separation in the overall gravitational lensing. 
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Parameter 


value 


Mi 


2.7 • l(r 7 M = 0.090 M ffi 


Mo 

_ 


2.4 • 1O" 4 M = 80 M e 


M = a/Mi Mo 


8.05 • 1O~ 6 M = 2.7M e 


X = Mo /Mi 


890 


^ it 


902 nas 


^e(M) 


6.5 nas 


#eff 


71.4 nas 


(X 


0.707 


a* 


2.983:° 


(tev)(«*) 


11.1 day 


(5m bs)(a*) 


0.083 



Table 2: Empirical values of microlensing parameters 

Our power law slope a* = 2.98 is steeper than the Salpeter value of 2.35. There is 
no reason to expect the two values to agree, since they are likely to be dominated by 
different processes in star and in primordial rogue planet formation. 

5. Summary and Conclusions 

With the quasar microlensing detection of a cosmologically significant fraction of 
the missing baryons (baryonic dark matter) now confirmed by high-cadence MACHO 
searches to the galactic center and to the LMC (Sumi et al, 2010, 2011), we investigate 
methods to permit investigation of the mass distribution function of the rogue planet 
population. The only theory that has predicted that the baryonic dark matter should 
exist as a population of rogue planets (Gibson, 1996) predicts that at time of formation 
their mass was 10~ 7 M , and that they should be found primarily in primordial Jeans 
clusters of approximately 6 • 1O 5 M (Nieuwenhuizen, Schild and Gibson, 2011). Such 
Jeans clusters have also been found in quasar milli-lensing at significant optical depth 
(Mao and Schneider, 1998). In the high-temperature and high-density universe at 
z ~ 1000, these sticky hydrogen spheres were formed by the usual void-condensation 
separation process and would have immediately interacted to form pairs, triples and 
pairs-of-pairs to eventually accumulate in an accretional cascade that quickly produced 
the first stars and left behind rogue planets of larger mass. Most Jeans clusters should 
have remained dark, but rapidly form stars when disturbed. A prediction of this theory 
is that all or most stars should be binaries. 

We also demonstrate how the expected microlensing signal apparent in observed 
brightness fluctuations in quasars can be analyzed with a statistical theory devised 
by RS91. We find that the observed amplitudes of the microlensing fluctuations are 
approximately a factor 2 smaller than predicted, and that this mimics the amplitudes 
expected if the particles are not all of the same mass, but instead are distributed 
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according to a mass function with slope approximately 2.98_ ; 5 , where the Salpeter 
slope is 2.35. 

Our study thus far has been applied to the Q0957+561 A,B quasar system (the 
First Lens) for which a great deal of information on the luminous quasar structure 
has been inferred in SLR06. However this radio source quasar in the Lo-Hard spectral 
state is probably not optimum since its luminous inner accretion disc edge is probably 
finer than the High Soft radio quiet objects, which are also more common. Thus in the 
future we will extend our study to the several radio-quiet lensed quasars with time delay 
measurements and significant observed microlensing residuals. 
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